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Summary 

The acoustic response of a semireverberant enclo- 
sure with two interacting, velocity-prescribed source 
distributions was analyzed using standard modal 
analysis techniques with a view toward a better un- 
derstanding of active noise control. Different source 
and enclosure dimensions, source separations, and 
single-wall admittances were studied over represen- 
tative frequency bandwidths of 10 Hz with source 
relative phase as a parameter. Results indicate that 
power radiated into the enclosure agrees qualitatively 
with the spatial average of the mean square pressure, 
even though the reverberant field is nondiffuse. De- 
creases in acoustic power can therefore be used to es- 
timate global noise reduction in a nondiffuse semire- 
verberant environment. As might be expected, para- 
metric studies indicate that maximum power reduc- 
tions of up to 25 dB can be achieved when secondary 
and primary sources are compact and closely spaced. 
Although less success is achieved with increasing fre- 
quency and source separation or size, significant sup- 
pression of up to 8 dB still occurs over the 1 to 2 Hz 
bandwidth. 

Introduction 

Concern about high levels of low-frequency inte- 
rior cabin noise for proposed advanced turboprop air- 
craft and recent developments in microprocessor con- 
trol technology has stimulated new interest in active 
suppression of aircraft interior cabin noise. If active 
suppression is effective for interior cabin noise, then 
bulky passive absorbers or heavy barrier-type ma- 
terials can be avoided, and weight can be reduced. 
The success of active suppression has been limited 
in three-dimensional acoustic spaces. Nevertheless, 
the modest success to date offers hope that substan- 
tial improvement may be obtained with currently 
available control technology through a better under- 
standing of interactions between acoustic sources in 
a semireverberant environment. To realize the po- 
tential improvement, additional analytical work is 
required which should provide an understanding of 
the phenomenon involved. An overview of previous 
experimental and analytical work is given in refer- 
ences 1 and 2. This earlier work is summarized in 
the paragraphs which follow. 

For the most part, past work on active suppres- 
sion has been concentrated on one-dimensional sys- 
tems. After Lueg demonstrated in 1934 that sup- 
pressions up to 25 dB can be achieved in a duct us- 
ing a monopole control source (ref. 1), other work- 
ers, notably Jessel and Swinbanks, investigated more 
complex control-source configurations such as dipoles 
and tripoles in an attempt to improve system stabil- 


ity and bandwidth. These efforts indicate that, in 
principle, significant active suppression over useful 
bandwidths can be achieved. With the recent emer- 
gence of microprocessor-controlled adaptive control 
systems, it appears that practical, active suppression 
for one-dimensional systems will be available in the 
present decade. 

The literature on active suppression for two- and 
three-dimensional acoustic spaces is sparse. In 1956, 
Conover (ref. 3) reported results of an attempt to 
suppress far-field radiation from a 15 000-kVA trans- 
former. He achieved suppression of 6 to 25 dB over 
a polar angle of about 23° and out to a radial dis- 
tance of 125 ft by mounting a single control source 
(loudspeaker) on the side of the transformer struc- 
ture. Beyond the 23°, no suppression was observed 
for regions where the original noise level was rela- 
tively high. In regions where the original noise level 
was low, the levels were enhanced. Another exper- 
iment on active suppression applied to transformer 
noise is reported in reference 4. In this experiment, 
noise reductions obtained at six positions inside an 
office (control space) located across a courtyard from 
a twin-transformer installation were measured. As in 
the Conover experiment, a single control source was 
located in the vicinity of the transformers. However, 
the experiment differed from the others in that build- 
ings either flanked or obstructed the direct transmis- 
sion path to the control space. At the transformer 
fundamental frequency of 100 Hz, no less than a 6- 
dB reduction was observed throughout the office; in 
some locations, up to a 25-dB reduction was ob- 
served. At the second and third harmonics, reduc- 
tions of up to 10 dB were localized to volumes of 3-f’t 
radii or less. 

The results reported by Conover (ref. 3) and 
Ross (ref. 4) are encouraging, because noise reduc- 
tions were achieved in three-dimensional acoustic 
spaces of practical interest in which offending sources 
were noncompact. Reductions decreased and became 
more localized with increasing frequency, as would be 
expected for a single control source in the presence of 
an extended or noncompact source. However, mul- 
tiple control sources, together with precision control 
of phase and amplitude may help compensate for de- 
creasing compactness with increasing frequency. 

Recently, an active control experiment was con- 
ducted by Zalas and Tichy (ref. 5) on a fully outfit- 
ted, modern turboprop aircraft. The main thrust of 
this experiment was to demonstrate that significant 
broadband control can be obtained in the vicinity 
of a control sensor located in an aircraft cabin. A 
single control source was used and the aircraft was 
operated at cruise conditions. Overall engine har- 
monic levels at the control microphone were reduced 


by up to about 23 dB at frequencies ranging from 
50 to 500 Hz. However, the volume extent of the re- 
duction was disappointing. Detailed intensity map- 
ping at the fuselage walls indicated the presence of 
acoustically noncompact sources of energy inflow and 
outflow which changed location with frequency. 

The experiment described in the preceding para- 
graph and other recent work suggest that success- 
ful active suppression in a three-dimensional acous- 
tic environment generated by noncompact acoustic 
sources requires judiciously located multiple control 
sources. Also, a real-time adaptive control system ca- 
pable of optimizing relative phasing and amplitudes 
using data from multiple control sensors is necessary. 
The acoustic interactions between the sources and 
environment are sufficiently complex that experimen- 
tal work cannot adequately explore all parameters of 
interest. Analytical guidance would therefore be of 
great use in ascertaining those parameters of greatest 
import ance. 

Using modal analysis, the interaction of two 
sources radiating into an absorptive enclosure was 
investigated. The parameters used were source and 
enclosure size, source separation and relative phase, 
and wall admittance. This study is an attempt to 
identify the relative importance of these parame- 
ters without addressing the transient behavior of the 
acoustic system. 
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arbit rary const ant s 

source width 
source half-height 
speed of sound 

eenter-to-centor source separation 
distance 

frequency 

resonant frequency of (f.rn.n) 
mode 

mode coupling coefficient 

= yTi 

nondimonsional propagat ion 
constant. k mM /k 

real and imaginary parts of 
propagation (“oust ant 

free-space wave number, uj/c 

propagation constant 


kx • ky , k% 

separation constants 

L X . Ly. L z 

enclosure dimensions 

L m . n 

integers 

AN/Af 

modal density 

P m.n 

pressure associated with (m.n) 
mode 

p(x. y. z. t) 

acoustic field pressure 

Rrri.n 

reflection factor associated with 
(rn, n) mode 

T(t) 

solution of separated differential 
equation for time 

t 

time 

11 

prescribed source velocity 

Prn.n 

particle velocity associated with 
(m. n) mode 

f/o 

primary-source velocity 
amplitude 

u(x , y, z.t) 

acoustic particle velocity 

hUrn.n 

power radiated in the (rn,n) 
mode 

Wt 

total power radiated by source 
configuration 

X(x),Y(y),Z(z) 

solut ions of separated differential 
equat ions for spatial dimensions 

x. y , z 

spatial coordinates 

a 

normalized source separation 
parameter 

V 

normalized admittance 


real and imaginary parts of 
normalized admittance 

V 

ratio of secondary to primary 
source velocity. <t>, s /<t> p 

V 2 

Laplaeian operator 

K. 

ratio of source width to enclosure 
width, a/ L z 

^ m.n 

constant defined in text 


enclosure width and height, to 
length ratios, respectively 

V 

ratio of source half-height to 
enclosure height , b/ L y 
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Superscripts: 

(*') 

(r) 

* 


nondimensional spatial coordi- 
nates, x/L x , y/L y , and z/L z , 
respectively 

density of air 

primary and secondary source 
velocity distributions 

phase difference 

reflection-factor phase angle 

angular frequency, kc 

denotes imaginary part 
denotes real part 
complex conjugate 


Substituting this expression into equation (1) gives 


X" Y" Z" 1 T" 

~x + v + ~z ~ 


= 0 


(3) 


Splitting off the time variable and equating the re- 
sulting equations to an arbitrary constant, — k} gives 


X" Y" Z" , 2 
X + Y + Z ~ 


(4) 


and 


\_T" 


-k 2 


(5) 


The time dependence is determined by a second- 
order, ordinary differential equation given by 


r" + ( kc) 2 T = 0 (6) 


A bar over a symbol denotes a spatial value. 
Angle brackets (( )) denote a time average. 

Analysis 

The rectangular geometry depicted in figure 1(a) 
was chosen as a basis for the analytical model. The 
model consisted of a three-dimensional rectangular 
enclosure with two velocity sources on one wall- a 
primary or "offending** source and a secondary or 
"control” source, radiating toward an absorptive soft 
wall of uniform admittance. The acoustic behavior 
of the enclosure was investigated by using standard 
modal theory as given for acoustical systems in ref- 
erence 6. The main physical quantity of interest was 
the power delivered to the absorbing wall as a func- 
tion of frequency. Source separation, amplitude, rel- 
ative phase, geometric dimensions, and wall admit- 
tance were treated as parameters. Also, the rela- 
tionship between the power developed by the sources 
and the spatial average of the time-averaged or mean- 
square pressure was investigated. 


with the general solution 

T ~ Aexp(ifcct) + B exp(-ffcct) (7) 

Since only steady-state solutions are of interest, one 
of the preceding arbitrary constants can be zero. In 
this study, the constant A is chosen to be zero, so 

T = B exp(-ikct) = Bcxp(-iut) (8) 

where the excitation frequency u) has been intro- 
duced. Splitting off the x-dependence in equation (4) 
gives 



X" 

X 


= -A ; 2 + /e 2 


(9) 


Therefore, the differential equation for the x- 
dependence is 

X" + k 2 X = 0 (10) 

for which the general solution is 


Solution of Wave Equation 


X — A\ exp(ik x x) + B\ exp (~ik x x) (11) 


The acoustic field inside the enclosure is satisfied 
by the homogeneous wave equation 


Similarly, a choice of separation constants fc 2 and k% 
for the y- and ^-dependencies provides the equations 


v 2 p 


1 ^ 1=0 

c 2 dt* 


( 1 ) 


Y" + k 2 y Y = 0 (12a) 


subject to boundary conditions imposed at the wall 
surfaces and the prescribed source velocity distribu- 
tion. A separation-of- variables solution is made with 
the assumption that the pressure can be written as 

p(x.y.z.t) = X(x) Y(y) Z(z) T(t) (2) 


Z" + k 2 y Z = 0 (12b) 

with solutions 

Y = A 2 cos (k y y) + # 2 si n(k y y) (13a) 

Z = A3 co$(k z z) + #3 sin (k z z) (13b) 
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Finally, tho throe separation constants are related 
such that 

k^ = k x + ky + A~? (14) 

The normal component of the acoustic particle ve- 
locity at the rigid wall boundaries must be zero: it 
follows via the linearized momentum equation that 


dp 

dp 

dp 


dp 

dy I 

y= 0 (, y 

V=L y ' dz 

z = 0 

dz 


When these boundary conditions are applied to the 
general solutions given by equation (13b). the eigen- 
functions for the y- and ^-directions are given by 


V = /t 2 cos(fc ?y t/) (16a) 


Z = vos(k z z) ( 16b) 

and the eigenvalues for k y and k z are given by 

m 7r 

J y 

717V 




k? = 




(17a) 

(17b) 


Equation (14) can now be solved for the remaining 
separation constant k x to give 


A N 47 r f^L x L y Lz 2 7r / {Lx,L y + LxL z + L y L z ) 

A / r: 3 2 r 2 

+ 4 ( Lx + L V + ( 21 ) 

8c 

This formula can be restated for geometrically similar 
enclosures as a function of the dimensionless param- 
eters kL x (alternatively kL y or kL z )\ 


AN 

a/ 


™(fcL x ) 3 + (Ml+/ V - — (^x) 2 

7T 2 


+ U+ M ,+/Z 2 ) ( ^ ) 


1 


( 22 ) 


where /i] = L y /L x and = L z /L x . Clearly, the 
modal density of such enclosures does not scale lin- 
early with kL x (kL y ,kL z ). This observation implies 
that enclosure size effects cannot be characterized 
solely on the basis of a dimensionless frequency pa- 
rameter. 


Modal Pressure Solution 

Combining the various component solutions given 
by equations (8), (11), (16a), and (16b) and ab- 
sorbing the arbitrary constants into C m ,nand R m ji * 
gives 


K 


rnji 



W 7V \ 1 

UJ 



where the nondimensional parameter K nun has been 
introduced. This parameter can be interpreted as a 
nondimensional propagation (‘oust ant for the (m.n) 
mode of the enclosure when viewed as a transmission 
duct . 

In subsequent sections of this report, K mjl is 
correlated with the eigenvalues K? m n for the cor- 
responding rigid wall enclosure, which can also be 
written in nondimensional form as 



> / iwvy \ ( nivz \ v 

Prnji = b m .nCOS ( *° S [° X P \' lh m,nkx) 

T R-m.n exp( iK 7n jikx)]ox p( tujf) (23) 

for tho pressure in the (m,n) mode. The component 
solutions have been assembled in such manner as to 
allow the enclosure to be interpreted as a transmis- 
sion duct. (See ref. (7).) The G m .n is a coupling coef- 
ficient that quantifies source coupling into the (m,n) 
mode of the enclosure. The R m ji i s a mode reflec- 
tion coefficient that specifies the relative amplitude 
of the reflected (m,n) mode traveling in the nega- 
tive ^-direction. The R m%n and the G m ,n coefficients 
must now be determined. 

The axial particle velocity follows from the x- 
direction momentum equation to give 


The modal resonant frequencies of the rigid wall 
enclosure an' given by 


Vn 


1 dPtnji 

ikpc dx 


( 24 ) 




nij) 


.) ' ^ l\ (\rti,n 


(20) for harmonic time dependence. Therefore, 


The modal density for the case with the rigid wall 
enclosure is given by equation (0.5.12) of reference 6 
as follows: 


V ni „ * m ,n „ cos(m7rr/) cos(nTrs) 

(><’ 

X [exp {ikL x Krn,n£) ~ ^?n,n *‘xp( A:Lj- /\ m ,n O] (25) 
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where the nondimensional variables f = x/L x , 77 — 
y/L y , and f = 2/L 2 have been introduced. Since 
the admittance at the absorbing wall is spatially 
uniform, there is no modal coupling upon reflection. 
(See ref. 6, section 9.3.) Therefore, the normalized 
admittance presented to each (m, n) mode is defined 
by 


0 = pc 


Vm.n 

Pm.n 


£=i 


(26) 


Thus, 


Um.n _ Km.n exp (iKm, n kLx) — Rm.n exp( — iKm, 
Pm,n pc cxp(iKm.n kL x ) + Rm.n exp( — iKyn, 


Solving for R m ^ n yields 


kL x ) 

kL x ) 


(27) 


Rm.n - ^ m ' W f exp{2iK m , n kL x ) 
&m,n “t P 


(28) 


t + a/c) < f < (t + (<* + 1 )k 


The secondary source velocity is given by 


(31b) 


= sin 


where 


i( f ~ \ (a + 1)K )] sin [h { n - \ + ")] ( 32a) 


i--)<^<(i+^ 


t - (a + l)/c) < c < (\~ aK 


(32b) 


Note that 2a is the edge-to-edge source separation 
measured in source widths, and d (fig. 1(b)) is the 
center-to-center source separation, where d and a are 
related by 

d = (2a + l)a (33) 


At the source plane (£ = 0), equation (25) becomes 


Equation (28) specifies the mode reflection coefficient 
in terms of the normalized termination admittance 
/?, the mode propagation constant K m , n , and the 
nondimensional frequency parameter kL x . 

Source Specification 

The final constant to be determined is the modal 
coupling coefficient G m ,n, which includes the effects 
of the source geometry and velocity distribution. In 
the present study, the source velocity distribution 
chosen was a rectangular, simply supported plate vi- 
brating in its fundamental mode. Results were ob- 
tained for symmetric primary and secondary source 
locations about the centerline of the y-z plane. (See 
fig. 1(b).) The source configuration velocity distri- 
bution is specified by 


Uiv-s) = Uo[$p{r).$) + rexp(t'4>)$ s (?7 : f)] (29) 

where T exp(z^) specifies the relative amplitude and 
phase between primary and secondary source veloc- 
ities. The primary source velocity distribution is 
given by 


7T 

1 M 


’ n 

1 

M 

_K 

v ~ 2 - aK ) J 

sin 

> 

("“2 

+ u )\ 
(3( 


where k = a/L z and v = b/L y . (See fig. 1(b).) The 
nondimensional coordinates 77 and f vary over the 
ranges given by 



< T) < 



(31a) 


Um.n = cos(m7rr?) cos(n7rf)(l - R m , n ) 

pc 


(34) 

This equation gives the acoustic axial particle veloc- 
ity at the source plane for the (m,n) mode. Such 
modes can be used to represent the source velocity 
distribution and thereby satisfy the boundary condi- 
tion imposed at the source plane. A summation over 
the mode indices m and n to M and N, respectively, 
yields 


1 M N 

U{r) ,f) - — ll Km,nGm,nCOs(rMrri) 
m ~ 0 n ~ 0 

X COS(n7T^)(l - R m ,n) (35) 

The G m .n can be extracted by using the eigenfunc- 
tion orthogonality condition. The result is 

G m ,„= ff% t ) 

x cos(m7rn) cos(ri7r^) dr) d < ; (36) 


where U{r],$) is given by equations (29), (30), 
and (32a), and the appropriate integration limits are 
given by equations (31) and (32). Also. 


A m ,n = j (m = n = 0) 


= - (m / 0, n = 0, or m = 0, n / 0) 


and 


Am.n = 1 (m/o, n/ 0 ) 
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Power Into Soft Wall 


where 


The expressions for the power radiated by the 
sources need to be developed at this point. However, 
it is useful to first introduce some notation changes 
related to the concept of inode cut -011 frequency. It 
is convenient to calculate acoustic quantities of inter- 
est. such as power, as a function of frequency. For 
the (m.n) mode, excluding the (0.0) mode, the prop- 
agation constant K rn n is imaginary for sufficiently 
low-excitation frequencies. (See eq. (18).) For such 
values of K m - n , the (m.n) mode is said to be cut- 
off. (See ref. 6.) These imaginary values of K m n are 
denoted by 


Rmji — ihrn.n 



2 -i (37) 


As excitation frequency increases. K m n passes zero 
to take on real values. The frequency at which 
K’m.n — d is called the cut -on frequency for the 
(w.n) mode. Above this frequency, the mode prop- 
agates as a wave-like disturbance; below this fre- 
quency. the mode decays exponentially. Real values 
of Kmji are denoted by 



All expressions for power and spatially averaged 
mean-square pressure (p 2 ) undergo a change in char- 
acter upon passing through the cut -on frequency for 
a given mode. Therefore, in the paragraphs that fol- 
low, expressions for power and (p 2 ) are presented for 
a given mode in both cut-off and cut-on states. 

For a duct terminated by its characteristic imped- 
ance (infinitely long or no reflection), the cut-off 
modes do not carry power. However, if the duct is 
not so terminated, then the possibility exists for net 
power to be carried by the cut-off modes, albeit small 
relative to the cut -on modes. Rv integrating the 
intensity (1/2 R v(pu*)) over the absorbing wall for 
the (tn.n) mode, the following result is obtained for 
the power delivered by the sources to the absorbing 
wall for the (w. ti) mode in its cut-off state: 




V- 

V* 


R m ,n |/<Vw| sill l Pnt ,n (39) 


and 


\Rrn.n\ — 


(*Sn - 0 {i) f + {fW)‘ 

\J (K$n + 0rf + ([l{r)y 

x exp (—2 K^ n kL x ) 


( 40 ) 


ifrm.n — arc tail 


— arctan 


is t 1 ) /•j(i) 

— fi ( r ) 

/id'O . fj(i) 
rnji i h 

T^T 


(41) 


Examination of the above equations reveals that the 
maximum power carried by a cut-off mode occurs 
precisely at the cut-off* frequency (i.e.. K m%n = 0). 
The [lower well below the cut-off* frequency is neg- 
ligible because of the dominance of the exponential 
term, especially if the distance between t he source 
plane and termination L x is significant.. 

For a cut-on mode the power delivered to the soft 
wall is given by 


Wn 


where 




(42) 


\Rm,n\ — 


(a - x r) ) + 

\J (xttn + fi r) ) 2 + 


(43) 


If \R m , n \ approaches unity for a cut-on mode (i.e., 
perfect reflection), there is no power transmitted no 
matter how large the coupling coefficient O m .n- On 
the other hand, if there is no reflection {R m ji = 
0), the duct carries power only in the positive x~ 
direction. (Generally, for 0 < R m ji < E there may 
be significant, interaction between the source and re- 
flected waves which can. in turn, affect v ^ a t-he 

term (1 — R m ji) m the denominator of the expres- 
sion for (See eq. (36).) Therefore, increasing 

the reflection coefficient may increase the net power 
flow into the enclosure for particular combinations ol 
source configuration, excitation frequency, and ab- 
sorbing wall admittance. 

Although acoustic power is a convenient mea- 
sure of overall acoustic coupling efficiency of a given 
source configuration relative to an acoustic environ- 
ment. it is the acoustic pressure field that gives rise 
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directly to the sensation of noise by the human ear. 
On the other hand, the acoustic pressure may vary 
dramatically with excitation frequency and receiver 
position in a scmireverberant environment at low 
modal densities. The goal of this study is more a 
qualitative understanding of source-enclosure inter- 
action than a detailed study of localized noise re- 
duction in such an environment. It is of interest, 
however, to compare a measure of spatially averaged 
sound pressure in the enclosure with power developed 
by the sources for both cut-off and cut-on modes as 
follows: 


(. P 2 ) = r r - T [ [ [ PP* drj d$ (44) 
L x L/yL z Jo Jq Jo 

For the (m, n) mode in its cut-off state, the result is 

(P 2 ) = ~\Gm. n \ 2 t 4- [«p(-aK$ w fcL,) - 1] 

4 -2 K&\ n kL x 

+ 2 R e(Rm.n) + - 1 - [exp( 2K&] n L x ) - 1] 

2 K^nkLx 

(45) 

where R m , n - \R m ,n\ exp(z is defined by equa- 
tions (40) and (41). For the (m. n) mode in its cut -011 
state, the result is 


(P 2 ) — ~\ Gm,n(* 
4 


2\Rm,n | sin(2A m,nkL x 4- 'ipm.n ) — sin (ipm.n ) 
2K ( Z ] n kL x 




(46) 


where \R m ,n\ is defined by equation (43) and rpm,n is 
given by 


m.n ~ arc tan 


-/?(»> 

- 0 {T) 


- arc tan 


K^] n +0 (r) . 


(47) 


Equation (46) is characterized by a middle term that 
approximates sin(2K m ^ n kL x ) / (2Km.nkL x ) , since 
is typically small. If the (m.n) mode is be- 
ing driven well above cut -on, the middle term will 
likely not be a significant contributor. 

In summary, expressions have been developed for 
calculating the spatial average of the mean-square 
acoustic pressure and the power delivered to an ab- 
sorbing enclosure for two symmetrically placed, rect- 
angular velocity sources. The results have been pre- 


sented in a form that emphasizes the dependence 
on the frequency parameters kL x . kL y . and kL z \ a 
frequency parameter involving the source separation 
distance kd: the normalized absorbing wall admit- 
tance ft: and the secondary or primary source velocity 
relative amplitude and phase T exp(z</>). Representa- 
tive results are presented in the following section. 

Results and Discussion 

It is convenient to use total power reduction as 
a measure of noise reduction in the semireverberant 
enclosure; however, the expression for spatially aver- 
aged mean-square pressure contained a term which 
precludes perfect correlation with the expression for 
power. Therefore, an effort was made to establish the 
degree to which correlation exists between reductions 
of power and spatially averaged mean-square pres- 
sure. This relation was studied at both the low and 
high ends of the frequency range. A positive corre- 
lation between power and spatially averaged mean- 
square pressure in a semireverberant enclosure can- 
not hold in the limit of a fully reverberant enclosure, 
because there would be no power radiated by the 
sources, but mean-square pressure would certainly 
not be zero. 

Explicit comparisons of power and spatially av- 
eraged mean-square pressure were calculated for en- 
closure dimensions of (L x . L y . L z ) — (10. 9. 8) in. 
source dimensions of (a, 2b) — (0.3, 0.6) m, and an 
absorbing wall admittance of O.lpc. The plots in fig- 
ures (2) and (3) show comparisons between power 
and spatially averaged mean-square pressure for 10- 
Hz bandwidths centered at 45 and 205 Hz. respec- 
tively. The modal densities for these two frequencies 
are approximately 0.8 and 11 modes per hertz, re- 
spectively. (See table I.) The top three curves of each 
figure labeled ”0°* : . “ref." and “180°" show the spa- 
tially averaged mean-square pressure (p 2 ) for three 
different source configurations consisting of the pri- 
mary and secondary sources driven in phase, the pri- 
mary source alone, and the primary and secondary 
sources driven 180° out of phase. Note that (p 2 ) 
is expressed in decibels referenced to the calculated 
value of (p 2 ) for the primary source alone driven at 
the initial frequency. The bottom three curves of 
each figure likewise show the power level developed 
by the same source configurations expressed in deci- 
bels referenced to 1 watt. 

The shaded areas in each figure show what might 
be achieved in terms of reduction in spatially aver- 
aged mean-square pressure or power reduction for a 
secondary source 180° out of phase with the primary 
source. Also, “turning on" the secondary source at 
zero phase causes both (p 2 ) and the power to in- 
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crease by about 6 dB for the low-frequency (short- 
wavelength) case, which would be expected for a dou- 
bling in area of a compact source. At higher frequen- 
cies (shorter wavelengths) the configuration exhibits 
characteristics of extended sources. The main pur- 
pose for these two figures, however, is to compare the 
spatially averaged mean-square pressure inside the 
enclosure with the power developed by the sources. 

A comparison of the power reduction and (p 2 ) re- 
duction can be ascertained by comparing the shaded 
areas between the reference and 180° curves of each 
figure', (dearly, there is excellent qualitative agree- 
ment for the 1()-Hz bandwidth centered at 45 Hz. For 
the 1()-Hz bandwidth centered at 210 Hz. the corre- 
lation is good, but there are some differences in the 
neighborhood of the sharp peaks. It will be shown 
in subsequent figures that the sharper peaks are as- 
sociated with cut-on frequencies of the "duct-like" 
modes of the enclosure. Also, a dramatic degrada- 
tion in the noise reduction is apparent in figure 3 
compared with that of figure 2. This degradation 
is discussed in more det ail subsequently. The focus 
of interest presently is that reductions in power out- 
put correlate well with reductions in (p 2 ) at both 
the low and high end of the frequency range. As 
discussed in the section "Analysis," the term involv- 

ing s\n(2K m j l kL x )/(2K„ } j l kL x ) in equation (40) is 
a potential contributor to (p 2 ) near mode cut-on fre- 
quencies. However, for the modal power, this term 
is absent, which may account for the localized differ- 
ences between (p 2 ) reduction and power reduction at 
the higher frequencies. To recapitulate, it appears 
that reductions in spatially averaged mean-square 
sound pressure' level correlate well with reductions in 
sound power level for modal densities from about 0.8 
toll modes pen* hertz in a semireverberant rectangu- 
lar enclosure. Therefore, the ('fleets of source separa- 
tion. size, phasing, and wall admittance obtained by 
studying only the changes in the sound power levels 
art' discussed in the sections that follow. 

Effect of Secondary Source Phase and Source 

Separation 

The effects of incremental phase changes of the 
secondary source on the combined power developed 
by several primary and secondary source configura- 
tions are examined in this section. Figures 4(a). to 
4(d) show the effects of secondary source phase 1 incre- 
ments of 30° on the combined power output of the 
primary and secondary source 1 configuration. Each 
plot shows power level in dB referenced to l W ver- 
sus frequency for secondary source phase parameter 
ranging from 0° to 180° for a given source 1 separa- 
tion. At the top e>f each figure* is shown the range of 


the separation parameter kd for the particular source 
separation shown in the sketch on the right-hand side 
of the figure. The figure includes the soft-wall ad- 
mittance value and the plot symbols for the phase 
parameter. The frequencies of the power response 
maxima correlate with the hardwall enclosure modal 
resonant frequencies, as listed in table II, that were 
calculated from equation (20). Included in the figures 
are the hardwall enclosure mode order designations 
associated with these maxima. 

In particular, figure 4(a) shows the effect of sec- 
ondary source phasing on total output power when 
the source center-to-center separation is set to the 
minimum of one source width for the smallest source 
width used in this study. (See table I.) For phase 
differences ranging from 0° through 150°, the power 
output is dominated by three maxima. These max- 
ima are labeled by modal indices corresponding to 
the hardwall enclosure modal resonant frequencies 
of 41.85. 42.94, and 46.25 Hz (see table II). As 
shown in figure 4(a), these frequencies correspond 
closely to the frequencies at which the maximum 
power responses of the semireverberant enclosure 
occur. For a secondary phase difference of 180°, 
however, the character of the power response curve 
changes dramatically, because the even-mode com- 
ponent associated with the 2 -axis direction cannot, 
be driven with an asymmetric source velocity dis- 
tribution. The correspondence between rigid enclo- 
sure modal resonances and the soft-wall enclosure 
modal behavior is expected, because t he admittance 
is chosen small to make the enclosure highly rever- 
berant. Also, the “break-even" phase angle (defined 
in discussion of fig. 5) at which the combined output 
power of the primary and secondary source config- 
urations is very nearly equal to the primary source 
power alone is about 120°. This phase angle arises 
from the effects of the interaction of two compact 
sources in close proximity, and appears mathemat- 
ically as 2h()<F p $, s cos 0 when the quantity UU* is 
calculated using equation (29). This quantity en- 
ters the power calculation through the factor |(» mji | 2 
in e(piat ion (39). When the secondary and primary 
phase differences are 180°, the minimum and maxi- 
mum power reductions across the entire 10- Hz band- 
width are about 12 and 34 dB. respectively. This 
substantial noise reduction is consistent with the 
small source separation relative to the wavelength 
(0.22 < kd < 0.27) and the source compactness 
(0.22 < ka < 0.27). 

Figure 4(b) shows the power developed by the 
sources at secondary source phase increments of 30 ( 
for a source separation of 9 source widths. 1 he 
key feature of this result compared with that of fig- 
ure 4(a) is the absence of a relatively uniform redue- 
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tion in power as the phase increases toward 150°. 
There is a substantial decrease in power between 41 
and 43 Hz. In the frequency range 43.4 to 49 Hz. 
the power generally increases with increasing phase 
difference. As shown in figure 4(a) these changes in 
the power with changing phase again correlate with 
the excitation of particular modes that can be iden- 
tified with hardwall enclosure modes. The (1.2.0) 
and (0.0.2) modes are symmetrical in the y— and 
z— directions and are therefore annihilated when the 
phase difference attains a value of 180°. On the other 
hand, the (0,2,1) and (1,2,1) modes are asymmetrical 
in the z-direction and are accentuated as the phase 
increases toward 180°. 

The source separation in figure 4(c) has been in- 
creased to 17 source widths. The power reduction 
achieved for this separation is less than that for 
the case with a separation of 9 source widths and 
is likewise associated only with the annihilation of 
the (1.2.0) mode as the phase increases toward 180°. 
Again, as in figure 4(b). there is significant power in- 
crease over the frequency range 43 to 50 Hz with in- 
creasing phase, (dearly, the (0,2,1) and (1,2,1) modes 
couple to this source separation better than for the 
9-source-width separation shown in figure 4(b). 

For the limiting source separation of 25.6 source 
widths (fig. 4(d)), the general trend is similar to 
those of figures 4(b) and (c). There is slight degra- 
dation in power reduction from 41 to 43 Hz, and the 
(1,0,2) mode is driven more strongly for this sepa- 
ration. Figures 4(a) to (d) summarize the detailed 
effects on total power output of dephasing the sec- 
ondary source with respect to the primary source in 
steps of 30° from 0° to 180°. In general, increas- 
ing source separation causes the enclosure mode re- 
sponse to change dramatically, which in most cases 
results in degradation of achievable power reduction. 
Figures 5(a) to (d) use the 0° and 180° plots from 
figures 4(a) to (d) to emphasize the output power 
reduction achievable relative to the power output of 
the primary source alone when the secondary source 
is driven 180° relative to the primary source. 

Sound Power Reductions for Secondary Phase 

of 0° and 180° 

Figures 5(a) to (d) show output power versus fre- 
quency for the same four source separations discussed 
previously, but each source configuration is operating 
in three different states, corresponding to the three 
curves shown in the figures. The top curve, denoted 
by the circles, corresponds to both sources operat- 
ing in phase; the curve denoted by the diamonds 
represents the power level for the primary source 
alone: finally, the curve denoted by the squares rep- 
resents the output power level for both sources op- 


erating 180° out of phase. For both primary and 
secondary sources “turned on.” the velocity ampli- 
tudes are equal, as indicated by the velocity ratio 
parameter T in the figure key. The curve represent- 
ing the power output for the primary source alone is 
taken as a reference or break-even power level, which 
the combined primary and secondary output must 
equal to obtain a payoff from an active control view- 
point. The shaded areas between the curves repre- 
sent the potential sound power reduction over the 
10-Hz bandwidth for the particular source separa- 
tion shown in the sketch, on the right-hand side of 
each figure for both 0° and 180° relative phase. A 
different type of shading corresponds to reductions 
achieved for each of the two control source phasings 
considered. 

Figure 5(a) shows that net power reductions rang- 
ing from 10 to 25 dB can be achieved over the en- 
tire 10-Hz bandwidth for the minimum separation of 
1 source width. In figure 5(b), the source separa- 
tion is 9 source widths. Clearly, a marked degrada- 
tion in the potential for sound power reduction oc- 
curs with this increased source separation. However, 
over a 1-Hz interval centered at 42 Hz. a reduction 
of up to 8 dB is still achievable. IJp to 3-dB reduc- 
tions also appear in the vicinity of 44 and 47 Hz. In 
figure 5(c), the separation is 17 source widths. For 
this separation, the frequency interval and magni- 
tude of power reduction centered at 42 Hz decreased; 
however, the regions of reduction associated with 
the two higher frequencies of 44 and 47 Hz in- 
creased and coalesced into a single continuous re- 
gion to yield power reduction ranging up to 8 dB. 
Finally, for the maximum separation of 25.6 source 
widths (fig. 5(d)). the frequency intervals of power re- 
duction become more limited with maximum reduc- 
tions of about 5 dB. Figures 5(a) to (d) demonstrate 
that small source separation gives consistent reduc- 
tions across at least 10-Hz bandwidths, whereas large 
source separation gives no systematic result across 
the same 10-Hz bandwidth. However, at frequen- 
cies where reductions do occur at large separations, 
the reductions may still be significant. Generally, 
as source separation increases, the magnitudes and 
bandwidths of the sound power reductions appear to 
decrease. 

Source Size Effects 

Effects of source size on power reduction were 
examined. Power level reductions for minimum and 
maximum source separations are presented for a 
10-Hz bandwidth centered at 45 Hz and a 10-Hz 
bandwidth centered at 205 Hz. 

Figures 6(a) to (c) show the effects of source size 
on power reduction for the secondary source driven 
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at a phase of 0° or 180°. The source separation 
parameter kd ranged from 1.98 to 2.47 for all three 
source sizes. Comparison of the shaded areas for 
figures 6(a) to (c) shows great similarity in the noise 
reductions as the source size changes, although there 
is some degradation in power reduction for the largest 
source size shown in figure 6(c). These results suggest 
that for values of kd between 2.0 and 2.5. power level 
reduction scales with kd. but is invariant with source 
size parameter ka up to a value of 2.5. 

Figures 7(a) to (c) show the effects of source size 
on power reduction for the secondary source driven 
at a, phase* of 0° or 180° over the frequency range 
of 200 to 210 Hz or over the kd range from 9.87 to 
10.87 and over the ka range from 1.10 to 10.22. No 
consist (mt reductions are observed in this frequency 
range for the three different source sizes. In fact., 
no significant power reduction is present except, for 
the minimum source separation (limited by edge 
contact) of figure 7(c). For the source separation of 
figure 7(c). both the separation parameter kd and the 
source* compactness parameter ka range from 9.87 
to 10.87. It appears that a modest power reduction 
can be obtained even for noncompact source's if the 
separation distance is minimized. However, for the 
large*!* sources the power reduction doe's not scale' 
with kd 

Effect e>f Enclosure Si/e 

Figure's 8(a) to (c) illust rate t he eficct of enclosure* 
size' on poweT mluct ion. Fe>r successive* demblings of 
the' enclosure eiime'iisiems, the range's of the* source 
size' parameter ka anel source* separation parame*te*r 
kd are* hole! cemstant although the encle>sure size 
parameter kL r assuine's a different range for each 
figure*. The* kL x range's are ele'picteel at the bottom e)f 
e*ach plot. The* minimum value of kL x is 7.8, anel the 
maximum value* is 88.0. As in the pre’ceeling figures, 
sketches of source* place'nu'iit relative te) the bounding 
('iiclosure* walls are ine*lueleel. For the source size* 
chosen. the* maximum separation of 5 semree wielths 
is de'termined by the smallest enclosure size. (Se*e 
liK- S(iv).) 

Comparison of the* power reduct iems for enelo- 
sures eli tiering in size* by a facte>r e)f four indicates no 
consistent tre'iuls ewer a source separatiem parameter 
kd range* e)f 1.1 centered on a value of appre>xiniately 
5.0. However, at those* frequencie's at which power 
dev reuses ele> e)e‘e*m‘. the*y range as high as 8 elB. 

Effect of Impedance Changes 

The rennaining figures de*al with (‘fleets of chang- 
ing the* absorbing wall admittance*. All the* previous 
results we're* base'd on a normalized admittance of 0.1, 


which corresponds to a relatively hard wall. Compar- 
ative results are presented for normalized admittance 
ft values e>f 0.1. 0.1 + O.li. and 0.2 for the interme- 
diate senirce size* (see table I) at the minimum and 
maximum senirce separat ion distance's. 

Figures 9(a) to (c) show the effect of increasing 
the normalized soft-wall admittance through the re- 
spective value's 0.1. 0.1 + O.h. and 0.2 at the min- 
imum source separation. Results are presented over 
the frequency range from 40 to 50 Hz. For the aehnit- 
tance increment, of 0.1 in the imaginary component., 
there are observable but relatively mine)]* (‘hange's 
in the shape's of the power re'duetion regions. (Sev 
fig. 9(b).) For the admittance increment e)f 0.1 in 
the real component (fig. 9(c)), there are significant 
changes in t he shapes of t he pe)wer mluetion regions, 
but a general similarity in behavior to figures 9(a) 
anel (b) is evielemt. Figure's 10(a) to (c) show the 
results for a similar computation for the maximum 
source separatum e>f 5 source widths. 

Again, there are relativc'ly minor changes in the* 
regions of power reduction fe>r the increment in the 
imaginary compement.. An increment of 0.1 in the 
real component , however, causes a substantial dc- 
e reuse in the power reduction. Thus, for tliis fn*- 
epiene'y range, a change of 0.1 in the real component 
of ft appears to have consistently greater effect than 
the same change in the imaginary component of ft. 

The remaining figures (11(a) to (c) and 12(a) 
to (<•)) show- the same sequence of parametric plots 
depicting effects of incremental change's in soft -wall 
admittance for the high end of the frequency range 
(a 10-Hz bandwidth centered on 205 Hz). For this 
frequency range, there is little change in the power 
reductions for a separation of 1 source width as 
the imaginary admittance component increases by 
0.1. (See figs. 11(a) and (b).) In figures 12(a) 
and (b). for the maximum source separation, there is 
again little change for the* increment in the imaginary 
admittance componemt. In contrast, figure 12(c) 
shows that, for the increment in the real component, 
there is no observable power reeluction achieved due* 
to tlu* secondary source for eitheT choice of phase. 
The results shown in figures 12(a) to (c) for the high 
frequency range thus exhibit trends similar to those 
seen in figures 11(a) to (c) for the low frequency 
range. 

Concluding Remarks 

Acoustic power radiated by two mutually inter- 
acting, rectangular velocity source's into a semirever- 
berant enclosure* was studied to determine' the rela- 
tive significance of controlling parameters as applied 
to active noise control. The range's of parameter val- 
ues were sufficient to conclude that the* power re- 
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duction is very sensitive to source size, separation, 
and phase. In general, the results obtained suggest 
that for low modal densities in a semireverbcrant 
enclosure, significant noise reductions are achievable 
for compact sources. However, for larger sources at 
higher modal densities, only modest power reduc- 
tions can be obtained, and these reductions degrade 
as the prirnary-to-secondary source separation pa- 
rameter increases. 

Source location symmetry considerations indicate 
that reduction in net radiated power to the enclosure 
is caused by changes in the radiation impedance on 
each source surface. These changes in mutual radi- 
ation impedance are directly related to changes in 
the acoustic pressure field in the enclosure and arc 
strongly dependent on source separation and phase. 
Reflections from the absorbing wall also affect the 
source radiation load; however, under semireverber- 
ant conditions, minor changes in absorbing wall ad- 
mittance are of secondary importance. 

These results suggest that active control of dis- 
crete frequency noise in a semireverbcrant enclosure 
depends on many interacting parameters. Power ra- 
diated by a primary or secondary source configura- 
tion varies markedly with changes in some parame- 
ters over critical ranges. Estimates of noise reduction 
achievable for aircraft cabin interiors requires calcu- 
lations for more realistic enclosure geometries and 
acoustic absorption characteristics. Also, source con- 
figurations involving the interaction of multiple pri- 


mary and control sources on opposite walls should be 
considered as a model for the acoustic transmission 
through aircraft fuselage sidewalls. 

NASA Langley Research Center 
Hampton. VA 23665-5225 
June 19. 1985 
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TABLE I. PARAMETER RANGES FOR FIGURES 4 THROUGH 12 




TABLE II. EIGENMODE FREQUENCIES FOR REVERBERANT 
ENCLOSURE ARRANGED IN MONOTONIC ORDER 

[L x = 10 m : Ly = 9 m : L z = 8 m : c — 343.5 rn / sec ] 


£ , m , n 

fi.m.n 

1,0,0 

" 17.18 

f 0, 1,0 

19.08 

f 0.0. 1 

21.47 

1.1.0 

25.67 

1,0,1 

27.49 

+0.1.1 

28.72 

1.1.1 

33.47 

2,0,0 

34.35 

+0.2,0 

38.17 

2. 1.0 

39.29 

2,0,1 

40.51 

1.2.0 

41.85 

+0.0.2 

42.94 

+0.2.1 

43.79 

2.1.1 

44.78 

1.0.2 

46.25 

+0.1.2 

46.99 

1,2,1 

47.04 

1. 1.2 

50.03 


^Zero-order, x-axis modes. 
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Figure 4. Continued. 
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(d) d = 25.6 source width: 
Figure 4. Concluded. 
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